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Abstract 

We determine the three-loop coefficient of the beta function in the 
asymmetric momentum subtraction scheme in Landau gauge. This scheme 
is convenient for lattice studies of a s , the running coupling constant of 
QCD. We present high statistics lattice results for a s in the SU(3) Yang- 
Mills theory without quark, compare with the three-loop running and ex- 
tract the value of the corresponding A_ parameter. We estimate the 
systematic error coming from four- loop terms. We obtain the result : 
A_ = 295(5)(15)^jfe^MeV. 
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Momentum subtraction schemes are interesting because they provide regular- 
ization independent renormalization methods which can be used non-perturbati- 
vely, contrary to the popular Ms schemes based on the dimensional regularization. 
Realized on the lattice, momentum subtraction schemes eliminate the need for 
lattice perturbation theory which has proved to be complicated and not very 
accurate. Usefulness of these schemes is exemplified in the non-perturbative 
renormalization of operators fl|] which becomes widely used by now in the lattice 
approach. 

In this paper we determine the three-loop running in an asymmetric momen- 
tum subtraction scheme, compare with our high statistics numerical results for 
a s and extract the value for the A_ parameter in the SU(3) Yang-Mills theory 
without quark. In section 1, we recall the definitions of the mom renormalization 

scheme and define g , a coupling constant convenient for lattice studies @, ||. 

In section 2, we give the relation between this coupling constant and some quan- 
tities computed in the ms scheme in ref.[[|. In section 3, we extract the three-loop 
coefficient of the beta function in the mom scheme and in section 4, we present 
our results for the coupling constant and the A_ parameter. 



1 Definition of the coupling constant in the mom 
scheme 

In this section, we suppose that we are able to compute or measure in some way 
the two- and three-point gluonic Green functions and define a convenient coupling 
constant. 

When one of the momenta is equal to zero, the Lorentz structure of the three- 
point Green function, G^ bc (p, —p, 0), can depend a priori on four tensors : 5p,uP P , 
fifipPv, bvpVp, an d Pp,PvVp- m Landau gauge, only one possibility is left by the 
transversality conditions and we can write : 

G$; bc (p,-P,0) = 2f^L u -^p p G^(p 2 ) (1) 

for the bare Green function. Note that only the completely antisymmetric group 
structure functions f abc can appear in this expression : the Bose symmetry of 
the three point function, i.e. the symmetry over p <-> — p, and a <-> b, 

combined with the oddness of the Lorentz tensor when p <-> — p forbids terms 
with the symmetric color structure d abc . 
The bare two-point Green function in Landau gauge writes : 

G$ ab ( P ,-p) = ^(<W-^r) G< 2 V), 
G$ ab (0,0) = 5 a % u G^(0). (2) 
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The scalars functions and can be extracted from the Green functions 
with : 

G (3 V) = f abc ^P p G^f c (P, ~P, 0) , 

G (2) (0) = t^E ^(0,0). (3) 
The wave function renormalization for the gluon, Z , and the renormalized 

MOM 

coupling constant in the mom scheme are then defined by: 



Z(^) = p 2 G^(p% 2 _ 2} (4) 



MOM 'P"=fi- 

,2> 



G( 3 V) Z 3 J^(fx 2 



2\ _ . MOM 



' "' G( 2 V)G( 2 V)G( 2 )(0) 



(5) 



where \i is the renormalization scale. The interpretation of eqs. ([§-[)]) is standard : 
the momentum scheme fixes the renormalization constants so that the two- and 
three-point functions take their tree values with the substitution of the bare 
coupling by the renormalized one. In the usual mom scheme, one chooses to work 
at the symmetric Euclidean point p\ = p\ = p\ = fj 2 . Here the mom scheme is 
defined from an asymmetric point, when one of the momenta is equal to zero. 



2 Connection with MS calculations 

In ref . [§]] , the two- and three-point Green functions have been computed at two- 
loop with the dimensional regularization in an arbitrary covariant gauge. They 
define the (bare) form factors from the (bare) vertex function: 

- Pms(<W 2 -^) Up 2 )}- (6) 

At zeroth order in perturbation theory, T\ = 1 and T2 = 0. 

If we restrict to Landau gauge and insert these form factors into the definition 

(H) for the coupling constant in the mom scheme, we get : 



MOM \ 47T / V 2 



Z*JL , (7) 
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where go is the bare coupling constant and a common scale, /x 2 , has been cho- 
sen for the renormalization point used to define g and for the dimensional 

MOM 

regularization. 

To renormalize in the Ms scheme, renormalization constants and a renormalized 
coupling constant have to be introduced. The definitions are: 

r /w3 _ P) o) = Z x I™" 3 (p, -p, 0) , (8) 
G^(p 2 ) = tfGVtf), (9) 

^0 - (^) £/2 ^^3 3/2 , (10) 

where the suffix Ms indicates renormalized quantities in the ms scheme. Z\ and 
Z 3 are the renormalization constants in the ms scheme; they are defined in the 
usual way to remove the 1 / e- 7 singular terms present in the bare quantities. From 
(HD we get Ti = Z{T\ and T 2 _ = Z X T 2 and from the definitions (§) and (§) : 

p 2 G^(p 2 ) = Z^p 2 G^(p 2 )\ = (11) 

Collecting everything, the relation between the mom coupling constant and Ms 
quantities is: 

where a = o 2 /47r and a = q^/Air. 

MOM MOM ^ MS 

3 The three-loop /? function in the mom scheme 

The f3 function in the mom and ms schemes are defined by : 

r\ ~ ~ ~ ~ 

a I \ mom A) 2 A 3 4 A 5 

p (a ) = \i — = at — -r-^or, —a; -at — ... 

MOM MOM Q^l 2% MOM 4:71 MOM 647H MOM 1287T MOM 

3(a) = u— ^ = -—a 2 -^a 3 - -^a 4 - (13) 

^MS V MS ' dfJ, 27T MS 47T 2 5e 647T 3 ^® V ' 

From 0, the perturbative expansions for the gluon propagator and the form 
factors 7]_ and T2 are known up to two loops in any covariant gauge. Eq.(|I 
can then be written as : 

fa_ «?_ 
1 + a + b —!£- + ..., I (14) 
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with known coefficient a and b. This gives a relation between the (3 functions : 
P^(a^) = /?_("-) 1 + a — + 3fo — ^ + ... . (15) 

MOM V MOM ^MS^ MS J \ 27T 167T 2 I 

Using eq.(|l4]) and the definitions ([13]), we can expand (|T5| ) in power of ct__ and 
identify the terms with the same power; we get : 

A) = Po, (16) 
Pi = Pi, (17) 
k = p 2 + 2{b-a 2 )p -4ap 1 . (18) 

Eqs. (|TE|,|r7D are usual results : the two first coefficients of the P function do 
not depend on the renormalization scheme. Eq. fllTf ) tells us that the two-loop 
results in Ms for the gluon propagator, T\ and T 2 (available in |||) are sufficient 
to get the difference between the three-loop coefficients of the P function in the 
mom and the ms schemes. In eq . fll8|) a, b and p 2 have to be computed in Landau 
gauge. But it is known that p 2 in ms does not depend on the gauge || so we can 
use the standard result in Feynman gauge, p 2 = 2857 - ^ N f + Nj @. 
Finally, the three-loop P function in the mcm scheme and Landau gauge is found 
to be: 

Po = 11 - |iV> , 
19 

Pi = 51- -N f , 
~ 186747 _ 1683 

Z35473 65 . \ Ar /829 16 . \ Ar2 16 AT , 

- {-it - T Ca ) w ' - (if " ^ Cs ) N > + t n " (19) 

where £3 ~ 1.2020569... is the value of the Riemann's zeta function and Nf is the 
number of flavor. In the following we will work in the flavorless case, Nf = 0, for 
which p 2 ~ 4824.31. 



4 A from a 

MS MOM 

The implementation of the mom scheme on the lattice is straightforward [^, |3|]. 
Varying the external momentum, the mom coupling constant can be obtained 
directly in one simulation for several values of the scale fi 2 . In ref.||, the two- 
and three-point Green functions have been measured on the lattice in Landau 
gauge with high statistics and for several lattice spacings and volumes in the 
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flavorless case. For large values of the momentum, lattice artifacts of 0(a 2 p 2 ) 
affect the Green functions. Numerically as expected these effects are seen to 
decrease when (3 increases. On the lattice the gauge fixing algorithm leads to the 
relation - sm(^ L )A fJ ,(p) = while p^A^p) does not vanish. And actually we 
have seen on our data [H| that the dominant part of 0(a 2 p 2 ) artifacts in a. is 

■- 1 MOM 

corrected by the substitution of the momenta p^ by the lattice momenta - sin(^| £i ) 
in eqs. ([!]) and @. This was already noticed in || ; in other contexts, authors 
have shown that this choice for the lattice momentum was indeed favored by their 
data, see for example f9], [7|]. 

We give now some results from the three-loop analysis, see || for details on the 
lattice settings. In Fig.ll], we give the behavior of the coupling constant, a. , as 

MOM 

a function of the scale /i obtained from the simulation and compared with the 
integration of the three-loop beta-function. A nice scaling is apparent for scales 
larger than 2 GeV. 

To calibrate the lattice at /3 = 6.2, we have used the value for the the lattice 
spacing which has been measured recently with a non-perturbatively improved 
action (free from 0(a) artifacts) 0. The measured value is : a _1 (/3 = 6.2) = 
2.75(18) GeV. Other lattices have been calibrated relatively to the one at j3 — 6.2 
with the results for a^fo ', the string tension in lattice units, published in M. We 



took: a-\(3 = 6.0) = 1.97 GeV, a" 1 ^ = 6.2) = 2.75 GeV and a~ l {(3 = 6.4) = 
3.66 GeV. 

Another way to exhibit the scaling is to extract a A parameter as a function of 
the scale and look for a plateau at high scale. Like in [|K]] we define in general 
the A parameter as : 



A = // exp 



-2vr 



f3 a — (// 2 ) 

MOM 



x exp < — 



MOM da 



MOM 1 

Ait 

1 2tt $! 



o 



(5{a) f3 ai (3 2 a 



(20) 



If we consider the expansion of the mom f3 function truncated at three-loop, (f20|) 
can be integrated to give : 



A 



-2tt 



uexp 

Mo: - J ' I (3 a (jj 2 



MOM 




x exp 



M - 4$ 



2/3 VA 



An 



1 + 



+ 



MOM 



arctan 



2f3 l +l3 2 a^/An 



2nf3 32vr 2 /3 
arctan 



(21) 



where A = 2/3 /3 2 — 4/3 2 (A > in our case). This can be consistently expanded 
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to give: 



0i 



A _ , „ exp | ^^_) x S ° fi + j o/ M 

\ P a {fi 2 )J \ 4vr J \ I6vr^ mom y 

(22) 

For each value of the scale u, there is a corresponding a and we can associate 

MOM 

an effective A parameter through the formulas above. This effective A pa- 

MOM ' MOM 

rameter should become a constant at sufficiently high fi 2 when the scaling settles. 
It should be noted that ( pif ) gives larger plateaux than (|2~2"1). But independently 
of the formula used to extract A, an estimate for the systematic error due to 
the influence of higher-order terms on this determination is needed. It can be 
obtained if we add by hand a four-loop term in the beta function and vary its 
coefficient (3% over a reasonable interval. Evaluation of these effects will be given 
below. 

The asymptotic ratio is known exactly p|: A_ = exp(— 1|) x A - 

~ 0.346 x A Let us forget for the moment the influence of higher order 

MOM 

terms and use the three- loop expression. In FigjJ we plot the effective Aqcd = 
0.346 x A as a function of u for several values of the volume and the lattice 

MOM 

spacing. Scaling manifests itself through the plateaux of Aqcd at large enough 
scale. The value of the plateau gives our "measurement" for the flavorless A_. 

From our largest physical volume ((3, V) = (6.0, 24 4 ) (for which V^ s ~ 2.4 fm) 

we obtain : A = (303 ± 5MeV). For comparison we give the results from the 

other lattices at (/3 = 6.0, V = 16 4 ), {(3 = 6.2, V = 24 4 ) and (p = 6.4, V = 32 4 ), 
all three with nearly the same physical volume, ~ 1.7 fm. We found A_ = 
314(3), 313(4) and 312(9) respectively. Scaling in (3 is striking and a comparison 
between the two different physical volumes shows that finite volume effects are 
moderate ( a few %). 

As explained previously, to estimate the systematic error from unknown higher 
order terms in the perturbative expansion we add a four-loop term in the beta- 
function and study the variation of A as a function of the unknown coefficient (3% 
varying up to some f3 l ^ m . We conventionally choose to limit this interval when 
the four-loop term in the beta function is equal in magnitude to the three-loop 
one at a = a hm = 0.4, namely (3 l ^ m = -7^/3 2 . We analyze the data and extract 
the values for A^ with 

MOM 

mom \ ft nfV (,,2\ / I 4tt I \ 16n/3^ MOM 




2/30/3x02 -m-Pifc , m -m 



2 



+ i — + 11 1 (23) 
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From the lattice with the largest physical volume this gives for our final result : 

A_ = 295(5)(15) a 1(/3 ~~ 6 " 0) MeV . (24) 
1.97GeV V ; 

where the first error comes from the statistics and the second one from the sys- 
tematics attached with the higher order terms. 
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Figure 1: The QCD coupling constant a as a function of the scale u (in 

MOM 

GeV). The full line is the three- loop running. 
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Figure 2: Aq C d (in GeV) as a function of the scale /i (in GeV) for different lattice 
spacings and volumes. 
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